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Abstract 

We study the Mott transition in Hubbard models with a degenerate band on dif- 
ferent 3-dimensional lattices. While for a non-degenerate band only the half-filled 
system may exhibit a Mott transition, with degeneracy there can be a transition 
for any integer filling. We analyze the filling dependence of the Mott transition and 
find that U c (the Hubbard interaction U at which the transition takes place) de- 
creases away from half-filling. In addition we can change the lattice-structure of the 
model. This allows us to study the influence of frustration on the Mott transition. 
We find that frustration increases U c , compared to bipartite systems. The results 
were obtained from fixed-node diffusion Monte Carlo calculations using trial func- 
tions which allow us to systematically vary the magnetic character of the system. 
To gain a qualitative understanding of the results, we have developed simple hop- 
ping arguments that help to rationalize the doping dependence and the influence of 
frustration on the Mott transition. Choosing the model parameters to describe the 
doped Fullerides, we can make contact with experiment and understand why some 
of the Fullerides are metals, while others, which according to density functional 
theory should also be metallic, actually are insulators. 
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1 Introduction 



The Hubbard model is an elementary model to study strongly interacting 
systems. In particular it can be used to understand the Mott metal-insulator 
transition. The original Hubbard model describes s electrons in a narrow band. 
With one electron per site, the half-filled system will become insulating for 
large enough correlation strength. To describe more realistic situations one 
has to generalize the Hubbard model to systems with degenerate (or near- 
degenerate) orbitals. Such degenerate orbitals can arise for example in molec- 
ular solids or transition metal compounds. For such a degenerate Hubbard 
model there will be a Mott transition not only at half-filling, but for all inte- 
ger fillings. It is then a natural question how the location of the Mott transition 
depends, for an otherwise unchanged Hamiltonian, on the filling. An other is- 
sue is the problem of how the Mott transition depends, for the same filling, on 
the lattice-structure of the system. To address these questions we have deter- 
mined the Mott transition for degenerate Hubbard models with various integer 
fillings and different lattice-structures. We have used Hamiltonians that de- 
scribe the alkali doped Fullerenes, since these materials have been synthesized 
in various integer dopings and crystal structures. 



2 Model and Method 



The inter-molecular interaction in solid Ceo is very weak. Therefore the energy 
levels of the molecule merely broaden into narrow, well separated bands [1]. 
The conduction band originates from the lowest unoccupied molecular orbital, 
the 3-fold degenerate t\ u orbital. To get a realistic, yet simple description of 
the electrons in the t\ u band, we use a Hubbard-like model that describes 
the interplay between the hopping of the electrons and their mutual Coulomb 
repulsion [2]: 

H = J]] ti m j m i C imcF Cj m i a + U nima-nim' a i . (1) 

(ij) mm 1 a i (ma)<(m' a') 



The sum (ij) is over nearest-neighbor sites. The hopping matrix elements 
Um,jm' between orbital m on molecule i and orbital w! on molecule j are 
obtained from a tight-binding parameterization [3,4]. The molecules are ori- 
entationally disordered [5], and the hopping integrals are chosen such that this 
orientational disorder is included [6]. The band- width for the infinite system 
is W = 0.63 eV. The on-site Coulomb interaction is U ~ 1.2 eV. The model 
neglects multiplet effects, but we remark that these tend to be counteracted 
by the Jahn- Teller effect, which is also not included in the model. 



2 



To identify the Mott transition we calculate the energy gap 



E g = E(N + 1) - 2 E(N) + E(N - 1), 



(2) 



where E(N) is the energy of a cluster of iV mo i molecules with N = n iV mol elec- 
trons (integer filling n). We determine these energies using fixed-node diffusion 
Monte Carlo [7]. Starting from a trial function \^t) we calculate 



where w is an estimate of the ground-state energy. The |tf( n )) are guaranteed 
to converge to the ground state |^ ) of H, if r is sufficiently small and \§>t) is 
not orthogonal to \$>o). Since we are dealing with Fermions, the Monte Carlo 
realization of the projection (3) suffers from the sign-problem. To avoid the 
exponential decay of the signal-to-noise ratio we use the fixed-node approxi- 
mation [7]. For lattice models this involves defining an effective Hamiltonian 
H c s by deleting from H all nondiagonal terms that would introduce a sign- 
flip. Thus, by construction, H cS is free of the sign-problem. To ensure that 
the ground-state energy of H c g is an upper bound of E , for each deleted 
hopping term an on-site energy is added in the diagonal of H e g. Since \^t) 
is used for importance sampling, H cf[ will depend on the trial function. Thus, 
in a fixed-node diffusion Monte Carlo calculation for a lattice Hamiltonian we 
choose a trial function and construct the corresponding effective Hamiltonian, 
for which the ground-state energy -Efndmc can then be determined without 
sign-problem by diffusion Monte Carlo. 

For the trial function we make the Gutzwiller Ansatz 



where the Gutzwiller factor reflects the Coulomb term U D = U "^ima^im'a' 
in the Hamiltonian (1). |$({7o)) is a Slater determinant that is constructed 
by solving the Hamiltonian in the Hartree-Fock approximation, replacing U 
by a variational parameter Uq. Increasing Uq will change the character of the 
trial function from paramagnetic to antiferromagnetic. This transition is also 
reflected in the variational energies obtained in quantum Monte Carlo, as 
shown in Fig. 1. Clearly, for small U the paramagnetic state is more favorable, 
while for large U the antiferromagnetic state gives a lower variational energy. 
Details on the character of the trial function and the optimization of the 
parameters can be found in Ref. [8]. 



= [1 - T {H - w)] n |* T ), 



(3) 



\^(U ,g))=g D mu )), 



(4) 
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Fig. 1. Dependence of variational (VMC) and fixed-node diffusion Monte Carlo 
(FN-DMC) on the trial function. Uo is the Hubbard interaction that was used for 
the Slater determinant in the Gutzwiller wavefunction ^t(R) = g D ^ $(Uo). The 
Gutzwiller parameter g has always been optimized. The results shown here are the 
energies (relative to the atomic limit) for a Hamiltonian that describes K3C60 (32 
sites), with U being varied from 1.25 (lowest curve) to 2.00 eV (highest curve). 

3 Results 

We now turn to the problem of the Mott transition in degenerate Hubbard 
models for different integer dopings and for different lattice-structures. The 
examples will be for integer-doped Fullerides A n C 60 , where A stands for an 
alkali metal like K, Rb, or Cs. Density functional calculations predict that all 
the doped Fullerides A n Ceo with n — 1 . . . 5 are metals [1]. Only Ceo and AgCeo 
are insulators with a completely empty/filled t\ u band. On the other hand, 
Hartree-Fock calculations for the Hamiltonian (1) predict a Mott transition 
already for U smaller than the band- width, and hardly any doping dependence. 
General arguments also suggest that the alkali doped Fullerenes should be 
Mott insulators, since the Coulomb repulsion U between two electrons on the 
same Ceo molecule (U ~ 1.2 eV) is substantially larger than the width of the 
ti u band (W 0.6 eV). It has therefore even been suggested that experimental 
samples of, say, the superconductor K 3 C 60 are metallic only because they are 
non-stoichio metric, i.e. that they actually are K 3 _5C 60 [9]. 

3.1 K 3 C 60 

In a first step we investigate what consequences the degeneracy of the t lu - 
band has for the Mott transition in K 3 C 60 . The analysis is motivated by the 
following simple argument [2,10]. In the limit of very large U we can estimate 
the energies needed to calculate the gap (2). For half-filling, all molecules will 
have three electrons in the t\ u orbital (Fig. 2 a). Hopping is strongly suppressed 
since it would increase the energy by U. Therefore, to leading order in t 2 /U, 
there will be no kinetic contribution to the total energy E(N). In contrast, 
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a) 



b) 



Fig. 2. Degeneracy argument: a) In the large-f/ limit hopping is suppressed for an 
integer-filled system, b) An additional charge can hop without extra Coulomb en- 
ergy. For degenerate orbitals and an antiferromagnetic background there are several 
different ways the extra charge can hop to a neighboring molecule. 

the systems with N ± 1 electrons have an extra electron/hole that can hop 
without additional cost in Coulomb energy. To estimate the kinetic energy 
we calculate the matrix element for the hopping of the extra charge against 
an antiferromagnetic background. Denoting the initial state with extra charge 
on molecule % by |1), we find that the second moment (l|i/ 2 |l) is given by 
the number of different possibilities for a next-neighbor hop times the single 
electron hopping matrix element t squared. By inserting J2j \ where \j) 
denotes the state with the extra charge hopped from site i to site j, we find 
(l\H\j) = V3t, since, with an antiferromagnetic background and degeneracy 
3, there are three different ways an extra charge can hop to a neighboring 
molecule (Fig. 2, b). Thus, due to the 3- fold degeneracy, the hopping matrix 
element is enhanced by a factor compared to the single electron hopping 
matrix element t. In the single-electron problem the kinetic energy is of the 
order of half the band-width W/2. The enhancement of the hopping matrix 
element in the many-body case therefore suggests that the kinetic energy for 
the extra charge is correspondingly enhanced. Inserting the energies into (2) 
we find that for the 3-fold degenerate system our simple argument predicts a 
gap 



E g = U-VSW, 



(5) 



instead of E g = U — W in the non-degenerate case. Extrapolating to inter- 
mediate U, it appears that the degeneracy shifts the Mott transition towards 
larger U. 

The above argument is, of course, not rigorous. First, it is not clear whether 
the result for E g that was obtained in the limit of large U can be extrapolated 
to intermediate U, where the Mott transition actually takes place. Also the 
analogy of the hopping in the many-body case with the hopping of a single 
electron is not rigorous, since the hopping of an extra charge against an an- 
tiferromagnetic background creates a string of nipped spins [11]. Nevertheless 
the argument suggests that orbital degeneracy should play an important role 
for the Mott transition. 
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Fig. 3. Finite-size corrected gap E g = E g —U/M—E g {U = 0) for increasing Coulomb 
interaction U as a function of 1/M, where M is the number of molecules. The 
calculations are for a Hubbard model with hopping matrix elements appropriate for 
K 3 C 60 . The band-width varies between W = 0.58 eV for M = 4 and W = 0.63 eV 
in the infinite-size limit. 

To check this proposition, we look at the results of the quantum Monte Carlo 
calculations for the model Hamiltonian (1) [2]. The Coulomb interaction U 
has been varied from U — ... 1.75 eV to study the opening of the gap. Since 
the Monte Carlo calculations are for finite systems, we have to extrapolate 
to infinite system size. To improve the extrapolation we correct for finite-size 
effects: First, there could be a gap E g (U = 0) already in the spectrum of the 
non-interacting system. Further, even for a metallic system of M molecules, 
there will be a finite-size contribution of U/M to the gap. It comes from the 
electrostatic energy of the extra charge, uniformly distributed over all sites. 
Both corrections vanish in the limit M — > oo, as they should. The finite-size 
corrected gap E g = E g — U/M — E g (U = 0) for systems with M — 4, 8, 
16, 32, and 64 molecules is shown in Fig. 3. We find that the gap opens for 
U between 1.50 eV and 1.75 eV. Since for the real system U — 1.2... 1.4 eV, 
K3C60 is thus close to a Mott transition, but still on the metallic side — even 
though U is considerably larger than the band-width W . This is in contrast 
to simpler theories that neglect orbital degeneracy. 



3.2 Doping dependence 



The degeneracy argument described above for K 3 C 60 can be generalized to 
integer fillings. Away from half-filling the enhancement of the hopping ma- 
trix elements for an extra electron is different from that for an extra hole. 
The effective enhancement for different fillings are given in Table 1. We find 
that the enhancement decreases as we move away from half-filling. Therefore 
we expect that away from half-filling correlations become more important, 
putting the system closer to the Mott transition, or maybe even pushing it 
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filling 


enhancement 


n = 3 


^3 « 1.73 


n = 2,4 


« 1.57 


n = 1, 5 


^ 2 +1 « 1.21 



Table 1 

Degeneracy enhancement for different integer fillings. 
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Fig. 4. Doping dependence of the Mott transition. The error bars indicate the 
estimate of the critical ratio U c /W for different integer fillings of the t\ u band. The 
calculations are for Fm3m Fullerides with fee lattice-structure and orientational 
disorder. The shaded region shows the range of U/W in which the doped Fullerenes 
are falling. 

across the transition, making it an insulator. We have analyzed the doping de- 
pendence of the Mott transition for the same Hamiltonian as used for K 3 Cgo, 
changing the filling of the t\ u band from n — 1 to 5. This model describes the 
Fm3m- Fullerides A n Ceo with fee lattice and orientational disorder [12]. The 
critical Coulomb interaction U c , at which, for the different integer fillings, the 
transition from a metal (for U < U c ) to an insulator (U > U c ) takes place, 
is shown in Fig. 4. As expected from the degeneracy argument, U c decreases 
away from n = 3. We note, however, that U c is asymmetric around half-filling. 
This asymmetry is not present in the simple degeneracy argument, where we 
implicitly assumed that the lattice is bipartite. In such a situation we have 
electron-hole symmetry, which implies symmetry around half-filling. For frus- 
trated lattices, like the fee lattice, electron-hole symmetry is broken, leading 
to an asymmetry in U c that is seen in Fig. 4. 



3.3 Dependence on lattice- structure 



To understand the effect of frustration in terms of the hopping arguments that 
we have made so far, we have to consider more than just one next-neighbor 
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Fig. 5. Basis states for four electrons on a triangle in the limit of large Hubbard 
interaction U. They are generated e.g. by hopping the extra charge clockwise around 
the triangle. After three hops the extra charge has returned to the original site. It 
takes, however, an extra lap to restore the spins to the original configuration. 

hop. The simplest system where we encounter frustration is a triangle with 
hopping matrix elements —t between neighboring sites. In the single-electron 
case we can form a bonding state with energy E m[n = —2t, but because 
of frustration we cannot form an anti-bonding state. Instead the maximum 
eigenenergy is -E max = t. Hence frustration leads to an asymmetric 'band' of 
width W = 3t. 

In the many-body case the situation is different. Like in the degeneracy ar- 
gument, we look at the hopping of an extra electron against a (frustrated) 
antiferromagnetic background in the large-C/ limit. For simplicity we assume 
a non-degenerate system, i.e. there is one electron per site on the triangle, plus 
the extra electron. In this case we have to move the extra charge twice around 
the triangle to come back to the many-body state we started from (cf. Fig. 5). 
Thus in the large-?/ limit the many-body problem is an eigenvalue problem of 
a 6 x 6 matrix with extreme eigenvalues ±2 1. In the degeneracy argument we 
have assumed that the kinetic energy of the extra charge is given by W/2. On 
the triangle, we find, however, that the hopping energy is by a factor of 4/3 
larger than that. This suggests that for frustrated systems the single electron 
band-width W in (5) should be multiplied by a prefactor larger than one. We 
therefore expect that frustration alone, already without degeneracy, shifts the 
Mott transition to larger U. 

To analyze the effect of frustration on the Mott transition we have determined 
the critical U for a hypothetical doped Fullerene A4C60 with body centered 
tetragonal (bet) structure, a lattice without frustration, having the same band- 
width (W = 0.6 eV) as the fcc-Fullerides, shown in Fig. 4. For U = 1.3 eV, we 
find a gap E g m 0.6 eV for the Fulleride with bet structure, while the frustrated 
fee compound still is metallic E g = 0. This difference is entirely due to the 
lattice-structure. Using realistic parameters for K4C60 [4] that crystallizes in 
a bet structure we find a Mott insulator with gap E g 0.7 eV, which is in 
line with experimental findings: E g = 0.5 ± 0.1 eV [13]. 
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4 Conclusion 



We have seen that, due to more efficient hopping, orbital degeneracy increases 
the critical U at which the Mott transition takes place. This puts the integer- 
doped Fullerenes close to a Mott transition. Whether they are on the metallic 
or insulating side depends on the filling of the band and the lattice-structure: 
Since the degeneracy enhancement is most efficient for a half-filled band, sys- 
tems doped away from half-filling tend to be more insulating. The effect of 
frustration, on the other hand, is to make the system more metallic. 
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